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Abstract. We present our photometric observations of GM Cep, a solar type vari- 
able in the young (~ 4 Myr) open cluster Trumpler 37. The star is known as a classical 
T Tauri star with a circumstellar disk and active accretion. GM Cep was suspected to 
undergo an outburst, thus a candidate for an EXor-type variable. In our monitoring 
campaign observations in 2010-2011, GMCep experienced a ~ 0.82 mag brightness 
decrease in the R band lasting for 39 days, and frequent, transient flare-like episodes 
with amplitude < 1 mag, each lasting for about 10 days. The brightening was ac- 
companied with a bluer color, presumably arising from increased accretion activity. 
Interestingly, the star also turned bluer in the fading phase. Combining the AAVSO 
and literature data, we found a quasi-cyclic peroid of ~ 3 11 days for the fading event. 
A possible mechanism for the fading could be obscuration by a clump of dust around 
the star We proposed that GM Cep therefore should be a UXor-type variable in the 
transition phase between grain coagulation and planetesimal formation process in the 
circumstellar disk. 



1. Introduction 

The Young Exoplanet Transit Initiative (YETI) is an international multi-site project 
which consists of a relay global network of telescopes (Fig. [T]l to search for exoplanet 
transit events in young open clusters (?). During an observing campaign, telescopes at 
different longitude range are used to monitor a young, nearby star cluster continuously 
for several days. Since 2009, the YETI has monitored two clusters, Trumpler (Tr) 37 
at an age of ~ 4 Myr (?) and 25 Ori at an age of ~ 10 Myr (?). While the primary 
goal of the YETI campaigns is to search for possible youngest exoplanets, the intensive 
monitoring produces data sets also valuable for young stellar variability study, which 
might be relevant to planet formation (?). 

The open cluster Tr37, at a heliodistance of 900 pc (?), is part of the Cepheus 
OB 2 association. Young stars are believed to evolve from Class II disk-bearing Clas- 
sical T Tauri Stars (CTTSs) to Class III diskless Weak-lined T Tauri Stars (WTTSs) in 
1 Myr to 10 Myr, when disk dissipation and planet formation are taking place. With a 
disk frequency of ~ 39% (?), Tr37 hence serves as a good target to search for and to 
characerize exoplanets in formation or early evolutionary stages (?). 
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Most pre-main sequences stars show irregular photometric variabihty. ? classified 
the variability into three categories. Type I variation is a periodic modulation caused by 
the rotation of a star with cool spots. Type II variation is caused by unsteady accretion 
or the rotation of hot spots on the star surface. Stars with Type III variation, also 
called UXors, with UXOri being the prototype, are hypothesized to suffer variable 
obscuration by circumstellar dust. 

GMCep (RA ^ 21 38 17.3, Dec ^ +57 31 23, J2000) is a solar type variable in 
Tr 37. It has a possible spectral type of G7 to KO, with an estimated mass of 2. 1 Mo and 
radius 3-6 RQ (?). The youth of GM Cep is exemplified by its emission-line spectra, 
prominent IR excess (?), and as an Chandra X-ray source (?), all characteristics of 
CTTSs. It has a circumstellar disk (?), with an accretion rate up to 10"^ Mq yr~\ 
which is 2-3 orders higher than the median value of the TTSs in Tr 37 (?). It is also one 
of the fastest rotators in the cluster, with v sin / ~ 43.2 km s"^ (?). 

Early observations of GM Cep by the YETI in 2009 suggested a possible EXor- 
type outburst. Previous studies on the light variability have been controversial. ? col- 
lected photometric data of GMCep from 1952 to 2007 in the literature, supplemented 
by their own intensive multi-wavelength observations. They suggested that GM Cep 
is an EXor type variable with an unstable disk and variable accretion rate. On the 
other hand, ? measured archival plates taken at Sonneberg and Harvard observatories 
between 1895 and 1993, and concluded the variability in the long-term light curve to 
be dominated by dips (possibly from extinction) superposed on quiescence states, in- 
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stead of outbursts caused by accretion flares. Here we present our multi-band optical 
observations of GM Cep, and discuss the possible mechanism underlying its variability. 



2. Observations and Data Reduction 

Our optical observations in the BVRc bands were carried out with the 1.0 m and 0.4 m 
telescopes at the Lulin Observatory (120.5E, 23. 3N) in Taiwan and the 0.81 m telescope 
at the Tenagra Observatory (110.5W, 31. 3N) in Arizona, USA. The Lulin One-meter 
Telescope (LOT) was equipped with a Princeton Instruments model 1300B (PI-1300B) 
camera, with a back-illuminated LEV CCD36-40 chip of 1340 x 1300 pixels. Each 
20 fim pixel corresponds to a plate scale of 0.51 "per pixel, yielding a field of view of 
ll'xll'. For the YETI campaign, a 0.55x focal reducer was used to enlarge the field. 
The 0.4 m SLT telescope at the Lulin Observatory was equipped with an Apogee Alta 
U9000 front-illuminated CCD camera with a Kodak KAF-9000 sensor. With 3056 x 
3056 pixels, each of 12 fj.m on a side, the field of view is 37' x 37'. Between 26 
August to 22 November in 2010, another camera was used instead. It was a back- 
illuminated Apogee Alta U42 with E2V CCD42-40, which has 2048 x 2048 13.5 /im 
pixels, giving a field of view of 28' x 28'. The 0.81 m Tenagra II telescope at the 
Tenagra Observatory was equipped with a back-illuminated SITe SI003 AP8p CCD, 
with 1024 X 1024 24 //m pixels. The field of view of the Tenagra II images is 15' x 15'. 
All the photometry images were reduced following standard routines to correct the bias, 
dark and flat-field. In addition, some data were taken with the 90/60 cm telescope of 
the University Observatory Jena (11.5E, 50.9N). For imaging, the telescope works in 
the Schmidt mode, with an effective diameter of 0.6 m. It was equipped with E2V 
CCD42-10 (STK), with 2048 x 2048 13.5 fim pixels. The images were subtracted by 
overscan and dark, then divided by the master flat. 

Aperture photometry was performed on GMCep by comparison with the seven 
reference stars within 3' of GM Cep used in ?. These comparison stars have comparable 
brightness as that of GM Cep. In each image we measured the FWHM of GM Cep and 
the comparison stars, then 3x, 5x, and 7x the largest FWHM among them were used 
for the aperture radius, inner sky radius, and outer sky radius, respectively. A linear 
regression between the instrument magnitudes and calibrated magnitudes of the com- 
parison stars was utilized to compute the brightness of GM Cep. Images with inferior 
sky conditions were excluded in the analysis. In addition to our own observations, vi- 
sual observations from the American Association of Variable Star Observers (AAVSO) 
from 2006 onwards were collected as well. 



3. Analysis 

The light curves of GMCep in the B, V, and R bands are depicted in Fig.[2]^a). Typical 
errors are smaller than the symbol sizes so are not shown. The variations in different 
bands follow well with each other, and abruptive fluctuations are obvious. Several days 
after our monitoring campaign started in 2010, GM Cep experienced a 0.82 mag fading 
in the R band lasting for about 39 days. The corresponding amplitude drop in the V band 
is smaller, about 0.68 mag. One also sees transient flaring-like episodes with amplitude 
less than 1 mag, each lasting for about 10 days. The nature of the variability is better 
diagnosed with color variations. Our data in the B band suffered large photometric 
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Figure 2. (a) The light curves of GM Cep in the B (open circles), V (gray circles), 
and R (filled circles) bands from our observations, (b) The R band light curve (filled 
circles) and V - R color curve (open circles) of GM Cep from our observations. 



errors, so the R band light curve and the V - R color curve are presented in Fig. |2jb). 
One sees that the color in general turned bluer as the star became brighter Interestingly, 
during the fading episode, GM Cep also had a bluer color 

We extended the time baseline of the light curve by including the data from the 
literature (?) and from the AAVSO archive. A quasi-cyclic fading period of ~ 318 days 
was discernible by eye inspection (Fig.|3]l. The NStED (NASA/IPAC/NExScI Star and 
Exoplanet Database) Periodogram Service was utilized to search for periodic signals 
with the Lomb-Scargle algorithm, after removing the long-term trend in the light curve. 
The first-ranked period is 311 days. 



4. Discussions 



The sporadic brightening episodes accompanied with a blue color can be accounted for 
by an enhanced accretion activity. The fading — a duration of 39 d with a blue color, 
and a period of 31 1 d — is puzzling. A possible mechanism could be due to obscuration 
by a clump of dust around the star, i.e., as suggested for UXor-type variables (?). We 
therefore suggest GM Cep to be an UXor, rather than an EXor (?). 

Assuming the clump is in Keplerian motion, and given the stellar mass of 2.1 Mo 
and an orbital period of P = 3 1 1 d, we derive the orbital distance ofa~ 1.16 AU. From 
the duration of the fading episode t, we have 



t 

P 



2Rc 
7.na ' 



(1) 



where f = 38.9 days, and Rq ~ 0.47 AU is the radius of the clump. The clump has a 
size of about 100 Ro. 

Furthermore, the amount of fading allows us to estimate the amount of obscur- 
ing dust along the line of sight, i.e., = 1.086 Nd Cd Qext, where Nd is the column 
density of the dust grains, crj is the geometric cross section of a grain, and 2^^? is the 
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Figure 3. The long term light curves of GM Cep in the B (open symbols), V (gray 
symbols), and R (filled symbols) bands. The data are from ? (squares), AAVSO 
(triangles), and our own observations (circles). 



dimensionless extinction efficiency factor. For particle sizes comparable to the wave- 
length, Qext ~ 1. Stars of a few Myr old would have the large grains settled into 
the midplane, so we assume the obscuration is attributed mostly to dust grains with 
an average radius of ~ 0.1 fim. Given the observed Ay = 0.68 mag, the column 
density of dust then is A^^ = 1.994 x 10^ cm"-^, which leads to a volume density of 
= Nd/2Rc = 1.424 x lO-^cm^l 
Since the clump is quite close to the star {a - 1.16 AU), we assume the composi- 
tion of the dust clump is mostly silicates, having an average density of p = 3.5 g cm~^. 
Therefore we derive the mass of the clump to be 

Md = ^7TR\dmd = 2.269 X 10^^ g = 2.269 X 10^^ kg, (2) 

which is about that of an asteroid. 

About a dozen UXors have been known thus far. Some show cyclic variability, 
with periods ranging from 8.2 days (?) to 1 1 .2 years (?). We propose that GM Cep be- 
longs to this category of objects experiencing highly variable circumstellar extinction. 
The mass we derived for the clump is only for the dust, and we have no evidence, even 
with a sufficient amount of associated gas, if the clump is on the verge of gravitational 
instability. Apparently a relatively small amount of dust mass along the line of sight to 
the circumstellar disk could cause the obscuration event we detected in our light curve. 
The intervening material does not need to shape Uke a clump, as long as the material 
is not evenly distributed azimuthally, e.g., a warped or arm-spiraled disk. Such inho- 
mogeneity in the young stellar disk might signpost the continuing process from grain 
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growth to the onset of planetesimal formation. Further characterization of the ununi- 
form disk of GM Cep, e.g., by polarization, reoccurrence of the bluing phenomenon in 
the fading epoch, infrared spectroscopy and submillimeter imaging in and out of the 
obscuration phase would shed more Ught on our hypothesis. 
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